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Abstract. Flow past a line vortex in a simple perfect fluid or superfluid gives rise to 
a transverse Magnus force that is given by the well known Joukowski lift formula. The 
problem of generalising this to multiconstituent superfluid models has been controversial 
since it was originally posed by the work of Iordanski in the context of the Landau 2- 
constituent model for ^He at finite temperature. The present work deals not just with 
this particular case but with the generic category of perfect multiconstituent models 
including the kind proposed for a mixture of Tie and ^e by Andreev and Bashkin. 
It is shown here (using a relativistic approach) that each constituent will provide a 
contribution proportional to the product of the corresponding momentum circulation 
integral with the associated asymptotic current density. 

1 Introduction 

For a simple perfect fluid with asympotically uniform density p say, the Magnus effect 
of a uniform background flow with relative velocity v % say in the rest frame of a vortex 
in the direction of a 3 dimensional unit vector t results in a force per unit length given 
by the well known (non-relativistic) Joukowsky formula as 

Ti = Kpe ijk £ J v k , (1) 

where k is the relevant velocity circulation integral. 

The question raised by Iordanskii |IJ of how this formula should be generalisd to the 
case of Landau's 2 constituent model for superfluid ^He at finite temperature has been a 
subject of controversy: the most widely accepted [| prescription is that of Sonin ||, |5|, 
but various alternatives have been proposed by other authors || [7|, ||. The present work 
clarifies the issue by demonstrating the existence of an elegant generalisation of the 
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Joukowsky formula ([!]) for an extensive class of perfect multiconstituent fluid models, 
including, as well as the Landau model, the Andreev Bashkin model || for a mixture 
of superfluid lie with ("normal") ^e. (However our analysis does not cover the more 
complicated subject of superfluid ^e). 

For technical convenience (not to mention the consideration that it is more accurate 
in contexts such as that of neutron star matter) the work is carried out using a (special) 
relativistic formulation . 



2 Perfect multiconstituent fluid dynamics 

As well as including the original Landau model in the (not so simple) Galilean limit Jl(J . 
two independent lines of early development of relativistic (and thus technically simpler) 
perfect multiconstituent fluid theory - using currents jQJ [12|] and momenta respec- 
tively as independent variables in (suitably constrained) variational formulations - were 
subsequently shown to be entirely equivalent ||14|| : the independent momentum covec- 
tors /i^ of the latter aproach are identifiable just as the dynamical conjugates of the 
independent conserved current vectors on which the former approach was based. In 
the latter approach, which is the most convenient for our present purpose, the funda- 
mental equation of state characterising a particular perfect multiconstituent fluid model 
is given by specifying the dependence of the relevant generalised pressure function \1/ 
on the independent momentum covectors /i*, and the associated current vectors 
are then obtained as the corresponding partial derivatives in the infinitesimal variation 
formula 

** = -XX*/**. (2) 

X 

For such a model, the complete set of equations of motion consists just of a set of 
particle and vorticity conservation laws of the form 

V v n£ = , n u x w ua = . (3) 

where, for any particular constituent with label x the corresponding generalised vorticity 
2-form is defined as the exterior derivative 

w* ff = V.vr* - , nl = //* + eM„ (4) 



of the generalised momentum covector 7r*, whose specification [T5|] allows for the possi- 
bility of coupling to an electromagnetic field 

F wr = V v A a -V a A v . V a F™ = 4irJ\ (5) 

with current and electromagnetic stress tensor given by 

J U = E e x < , T F \ = ^(F^F^ - \F^F, p gl) , (6) 
where e x is the electric charge, if any, per particle of the x th species. 
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In such a model a (flat space) conservation law of the usual form 

v„:r ff = o, 



(7) 



will be satisfied by the relevant total stress energy tensor, which takes the form 

T l = E n >l + W v + Tf v = + - J°A V + Tp u . (8) 

X X 

The transport law (|3p for the vorticities is such that they will remain zero if they 
are zero initially. Thus we shall have 

<, = 0, (9) 

not just for cases of superfluidity or superconductivity (i.e. cases for which the mo- 
mentum covector is the gradient of a condensate phase scalar) but even for "normal" 
constituents in configurations of the kind to be considered here, in which a perturbing 
vortex moves through an asymptotically uniform medium characterised by vanishing of 
the asymptotic background value (indicated here by an overhead bar) not just of the 
current (as is necessary for uniformity) but also of the electromagnetic field, and (in an 
appropriate gauge) of its vector potential, 

J u = F ua = 0, A a = 0. (10) 

This must necessarily be the case (the Meissner effect), with the implication that the 
uniform background value of the stress energy density tensor will be given simply by 

X 

whenever even just a single one of the uniform background constituents is supercon- 
ducting (since e x F ua = V„7r* - V a ix* - w^ a ). 



3 Specification of lift force on vortex 

The subject of this investigation is an asymptotically uniform vortex configuration that 
is stationary with respect to a rest frame characterised by a uniform timelike unit 
symmetry generating vector field k^ say, and that is aligned in the direction of a uni- 
form orthogonal spacelike unit symmetry generating vector field in a flat background 
spacetime using Minkowski coordinates. Provided the corresponding conditions of sta- 
tionarity and longitudinal symmetry apply to F va not just outside but even within the 
vortex core region, they will be applicable to the potential A a in a suitable gauge, and 
hence also, not just to the gauge independent covectors [i* but to the corresponding 
generalised momenta tt* as well. In view of the vanishing (§) of the vorticity vectors 
(f|), the condition that the momentum covectors should be invariant with respect to 
the action of the the uniform symmetry generating vector fields k u and l v can be seen 
to imply the uniformity of corresponding sets of generalised Bernouilli constants, 

k u 7T x u = k u W x u , rvr* = £ u Wl . (12) 
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The force per unit length, T v acting on such a stationary longitudinally invariant 
vortex can be evaluated as the integral round a circuit s say surrounding the vortex in 
an orthogonal 2-plane in the form 

Fu = <j> fvds , (13) 

where ds is the proper distance element given by ds 2 = g vu dx v dx a and /„ is the local 
force density that is given by 

f v = v a T° v , (14) 

in terms of the unit normal covector v a which will be given in terms of the antisymmetric 
background measure tensor e\^ U(T by 



4 Generalised Joukowski theorem 

It is to be observed that, as a consequence of the conservation law (0), it makes no differ- 
ence what circuit is employed for evaluating T v . We are thus allowed to choose a circuit 
sufficiently far out for reliability of our smoothed fluid description (whose physical va- 
lidity might be questionable near the core) to be ensured, and also for the deviation 
from the uniform background value T a v to be evaluated as a linear perturbation: 

T a u - T% = ST a u + 0{5 2 } . (16) 

Since the force integral for the unperturbed uniform background must evidently vanish, 
T v = by symmetry, the corresponding value in the presence of the vortex will be 
given by 

F v = 5F v + 0{5 2 } . (17) 

Using (D) and (JlTf) it can immediately be seen that the required first order variation 
will be given by 

5T° V = ]T (n<K + K^i - 9> P JK) ■ ( 18 ) 

X 

Using the decomposition of the 4-dimensional spacetime metric in the form g v a = 
rj l/ a + as the sum of the (rank 2) operators of projection respectively parallel to 
and orthogonal to the vortex given by rf a = —k u k a + £ u £ a , and -L^= % M! ^% M(J and us- 
ing the possibility of taking the Bernouilli constants (produced by parallel projection) 
outside the integration, (|18D provides a result expressible simply as 

8T V = (K * £ ^ 6CX + ^ 5 AJ , (19) 

X 

where for each species x the corresponding momentum circulation integral C x and 
current outflux integral _D X are defined by 

C x = l^ v dx\ D x = (fn°v a ds. (20) 
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The irrotationality condition (|9]) ensures that C x is independent of the choice of circuit, 
and the current conservation law ([5]) ensures that the same will apply to D x , which 
furthermore will simply vanish, D x = 0, provided there is no current creation in the 
vortex core. Thus by the fact that the uniform background value of the circulation 
integrals must also vanish, C x = 0, and by taking the limit in which the circuit is taken 
to a very large distance outside, one obtains an exact net force formula of the simple 
form 

^=^E C X- ( 21 ) 

X 

This result is the required (relativistic, multiconstituent) generalisation of Joukowsky's 
well known formula ([I]) for the single constituent case. What it means is that each 
constituent contributes an amount proportional to, but orthogonal to, its asymptotic 
current vector, with a coefficient given by the corresponding momentum circulation 
integral. 

5 Application to the Landau model 

The particular example that motivated this work is that of superffuid 4 He at finite 
temperature, as described by the Landau model in terms of just two constituents with 
conserved 3-dimensional current densities n l = n v l and n l = n„v l of which the first 

a a a f3 0/3 

represents Helium atoms, i.e. "dressed" alpha particles, characterised by a "rest mass" 
m a , and the second represents units of entropy, characterised by a vanishing rest mass 
nip =0. As in the less specialised case of the 2 constituent (zero temperature) limit of 
the Andreev Bashkin model for which the second constituent is 3 He with non vanishing 
rest mass, ~ 3m a /4, the Newtonian limit description can be formulated in terms of 
a total mass density and 3 dimensional mass current 

P = P a +P?, P i = py a + P p vi , (22) 

where p a = m a n a and p p = m^n^, so that the latter vanishes in the particular case 
of the Landau model. The total mass current is identifiable with the total momentum 
density pi = n a p1 + n^pl, in which, due to the effect of "entrainment" (which is 
describable in terms of "effective masses" that are different from the bare masses) the 
vanishing of the second contribution to the mass current does not imply absence of the 
second momentum contribution proportional to p\. 

The pseudo-velocity that is commonly referred to as the "superffuid velocity" is 
defineable by v s i = m" 1 /!". In the Landau case (unlike the generic Andreev Bashkin 
case) it is not possible to define an analogous pseudo velocity for the other consituent, 
because of the vanishing of , and the quantity commonly denoted as and known as 
the "normal" velocity is simply identifiable with the velocity of the entropy current, i.e. 
v n = v l~ ^ n a mass an d momentum decomposition of the commonly used (effectively 
"mongrel") form 

P = Ps + Pn , P l = Ps v s + PnK ' ( 23 ) 
the coefficients p s and p N must not be confused with p a and p (of which the latter is 
zero in the Landau case characterised by by p = p a ). In the Landau CcLSC, clS well as in 
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the generic Andreev Bashkin case, there are two independently conserved momentum 
circulation integrals, 

C*= (f ni dx i , C? = I p\ dx i , (24) 



of which, by the superfluidity property, the former is quantised, C a = h. Since m a ^ 
we can write 

C a = m a K s , C = rripKg + <j> — (v N i — v si )dx l , n s =j>v si dx\ (25) 

in terms of the pseudo- velocity circulation integral k s , which has no "normal" analogue 
in the Landau case, because of the vanishing of m* 3 . Thus (in the rest frame of the 
vortex) using the notation 'fey = Eijk£ k , one finally obtains a non-relativistic force 
formula of the form 

Ti = %j(K a p a vl + Cm^vp = K s p s %jvi + T xi , (26) 

where in this last version the first term is what is commonly referred to as the "superfluid 
Magnus force" contribution, while the remaining "Iordanskii" correction term is found 
to be given by 

^ = (C^ + /?s (p N -^))%^. (27) 

The third term in this expression is needed for the generic Andreev Bashkin case, but 
drops out for the special Landau case characterised by p p = 0. 

If, as well as setting p f3 to zero, one adopts the plausible supposition that the "nor- 
mal" circulation will vanish, C 13 = 0, then the first term also drops out so that our 
formula will reduce to a form that is in exact agreement with the result that was de- 
rived by Sonin|4], |5j and confirmed, on the basis of a more rigorous microscopic analysis 
of phonon dynamics, by Stone 0. 

This widely accepted conclusion has however been vigorously contested by Thouless 
and coworkers J7^ || who have used a more sophisticated - though not obviously more 
reliable - kind of microscopic analysis to argue that the Iordanskii force term T x % van- 



ishes, leaving just the purely "superfluid" term (namely fe^- K s p s v^) in (|26|) . As prima 
facie evidence in favour of this dissident conclusion, it is to be observed that in the 
limit when there is no relative flow at all (i.e. v* = = 0) then - as a requirement for 
compatibility with strict stationarity - the long term effect of the small "normal" vis- 
cosity contribution that was neglected in the preceeding analysis will impose a rigidity 
condition to the effect that v 1 = throughout. This imperative entails small deviations 
from strict irrotationality of the normal constituent except in the incompressible case 



for which the ratio p n /n fj is exactly uniform, and it ensures in any case by fl25D that 
the normal momentum circulation round a circuit at large distance will be given by the 
formula C 13 = (m g — p N /n )n s whose substitution in fl27D does indeed give, T x i = 0. 
However this simple counter argument is inconclusive because - as shown every time an 
ordinary light aircraft takes off - the stationary circulation value due to the long term 
effect of slight deviations from strictly inviscid behaviour will change as a function of 
the relative flow velocity. 
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To sum up, the present work shows how the Iordanskii force is given simply as a 
function of the (in the inviscid limit conserved) "normal" momentum circulation integral 
C 13 , but the issue of the appropriate value for this parameter in a realistic steady flow 
confguration is beyond the scope of a perfectly conducting fluid treatment such as is 
provided here. Experience with the analogous aerofoil problem in the context of aircraft 
engineering suggests that the final resolution of this issue may involve subtleties that 
have have eluded even the most sophisticated analysis available so far. 

We wish to thank Uwe Fischer, Michael Stone, and Grigori Volovik for instructive 
conversations. 
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